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Abstract. - We study the effect of mesoscopic fluctuations on a magnetic impurity coupled to a spatially 
confined electron gas with a temperature in the mesoscopic range (i.e. between the mean level spacing 
A and the Thouless energy iS-rii)- Comparing "poor-man's scaling" with exact Quantum Monte Carlo, 
we find that for temperatures larger than the Kondo temperature, many aspects of the fluctuations can be 
captured by the perturbative technique. Using this technique in conjunction with semi-classical approxima- 
tions, we are able to calculate the mesoscopic fluctuations for a wide variety of systems. For temperatures 
smaller than the Kondo temperature, we find large fluctuations and deviations from the universal behavior. 



The nature of many-body interaction effects at the nanoscale has been of great interest recently. 
Examples include, for instance, charging and correlation in quantum dots [1], superconductivity or 
ferromagnetism in metallic grains [2], and transport through effectively one-dimensional systems [3]. 
Since interference effects are ubiquitous at the nanoscale, it is natural to examine the interplay between 
interference and many-body interactions in these systems. 

A magnetic impurity coupled to an electron gas is a classic many-body problem, known as the 
Kondo problem [4]. In the simplest case, one takes a spin- 1/2 impurity coupled antiferromagnetically 
with an exchange constant J to an electron gas described by a constant density of states po within a 
bandwidth D. The physics has two distinct temperature regimes - at high T the electrons scatter off the 
impurity in inelastic spin-flip processes, while at low T there is only elastic scattering as the impurity is 
screened by electrons forming a bound state singlet. The crossover takes place at an energy scale T^, 
which in the weak coupling limit ~ De~ x l Jpa . All impurity properties (i.e. impurity contribution to 
thermodynamic quantities and impurity correlators) have temperature dependencies that are universal 
after rescaling with T^. This simplest case corresponds to a bulk piece of metal in which there is no 
interference. 

More recently, the effect of interference on Kondo impurities has been studied in various open 
systems. In disordered metals, the presence of magnetic impurities affects the weak-localization prop- 
erties [5,6], and the fact that each impurity sees a different local density of states may give rise to 
non-Fermi-liquid signatures in the metal-insulator transition [7]. In quantum point contacts, Friedel 
oscillations associated with the proximity of a boundary lead to a coupling between the magnetic im- 
purity and the electron gas which depends on both the position of the impurity and the energy of the 
electron state; as a consequence fluctuations of the Kondo properties are observed [8]. 

For closed systems, on the other hand, it is natural to ask: What is the effect of the confinement on 
the many-body physics? In ultra-small metal particles [2], for instance, the mean single-particle level 
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spacing A can be comparable to the bulk Kondo temperature T^. Clearly, a magnetic impurity in such 
a particle will show new physical properties associated with the new scale A produced by confinement. 
This has been studied theoretically in both metallic [9] and semiconductor heterostructure [10, 11] 
systems. However, most of these studies have neglected the fluctuations that are inherent in any 
mesoscopic sample (see however [12, 13]). Indeed, if A is the microscopic scale induced by the 
confinement, this latter brings in yet another energy scale, the Thouless energy £rh —h/t{ where tf 
is the time to cross the nanostructure. The range of energy between A and E^h is what is referred to 
as the mesoscopic regime. For a magnetic impurity coupled to electrons confined in a fully coherent 
nanoparticle at temperatures within this mesoscopic range, both the energies and the wave functions 
will display mesoscopic fluctuations which are stronger than those in open systems. 

The goal of this study is to understand the effect of these mesoscopic fluctuations on the bulk 
universal behavior. Experimentally, our work bears on two situations: first, an ultra-small metallic 
particle [2] containing a single magnetic impurity, and, second, the tunable Kondo impurities that have 
been formed using quantum dots [14, 15]. We return in particular to the latter at the end of the paper. 
Our approach is to combine analytic renormalization analysis [8] with exact numerical results. The 
combination allows a quantitative understanding for T> T^. 

To describe a magnetic impurity coupled to a band of confined electrons, we start with the single- 
particle electron states e a and wave functions <fi a (r) in the absence of the impurity. The impurity is 
modeled as a localized state at r = with onsite repulsion U. Finally, there is hybridization between 
the impurity state Cd and electron states c a ; for a <5-like interaction, the electron gas is completely 
described by the local density of states (LDOS) pi oc ( e ) = 2~Z Q 10a (0) | 2 <5(e — e a ). 

The combination of electron states that is most localized on the impurity is ?/v (0) = J2 a § a (®) c aa- 
The impurity couples to only this state of the electron gas. With locality and SU(2) symmetry, the 
simplest model is a symmetric (i.e. Cd = —U /2) Anderson model, 

H A nd = ffband + #imp + V^Ki^Lcd* + h.C.] (1) 

acr 

with i? band = J2 a ,a e aC f aa c aa and H imp = J2 a ^ d a c do + Und^ndi. We always take U large enough 
so that (1) is related through a Schrieffer- Wolff transformation to the Kondo model 

-ffKondo = -ffband + ^ 4^ - ^ a a x ~°*<Ti oi Cjiai (2) 
a/3 

where J a p = J(f)* a (0)^,3 (0) and J = 8V 2 /U. The band in the Kondo model must be cut at an energy 
scale D cut . (D cut = 0.18 U for the symmetric Anderson model [16,17].) Numerical calculations 
will be performed with the Anderson Hamiltonian Eq. Q, but most of our discussion will involve the 
Kondo form (0. 

Our goal here is to understand the fluctuations in measured quantities as the realization of the 
4> a (0) and e a changes. We focus on the case of chaotic confinement, for which a random matrix theory 
(RMT) description of the one-body physics is valid. For definiteness, we study the local susceptibility, 

X= [ dr(S z d (r)S z d (0)) (3) 
Jo 

where S z d (r) = e TH ^{rid\ — nd[)e~ TH . Experimentally, while this quantity is inaccessible for mag- 
netic impurities, we indicate in our concluding discussion how to measure \ directly using quantum 
dots. In addition, most of our analysis and conclusions can be extended to other properties. In the ab- 
sence of mesoscopic fluctuations (e a equally spaced and <p a ~ const.) and in the regime A « T C [/, 
X shows universal behavior [4]: T\ = /(T/T^-) so that the energy scale is the only parameter. 
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Fig. 1 - (a) QMC data for the local susceptibility of the Anderson model for the bulk-like case (open symbols) 
and three realizations of a mesoscopic electron sea (filled symbols). Also shown are the corresponding universal 
curves \ = I{T/Tk)/T for a realization dependent Tk obtained from Eq. (the correspondence between 
symbols and lines is the natural one). Inset: for realisations for which a good fit with a universal curve is obtained, 
comparison of the tJJRG extracted from this fit , with the T^ CIt from perturbation theory [Eq. Q]. (b) The 
smoothed local density of states at the impurity site for the same three realizations and the bulk case [line style 
matches the corresponding universal curves in (a)]. Note the striking non-universal curves for certain realizations 
and the correlation with variations in the density of states. The parameters are: V = nV 2 /A = 125A, <7 = 250A. 
T is scaled to the corresponding bulk Kondo temperature T% — 10. 2A (value from Bethe ansatz). The realization 
dependent (f> a and e a are drawn from RMT-GOE and for the bulk case are cj) a = 1 with e a equally spaced. 



Numerics. - Before turning to an analysis of the Hamiltonian Eq. (0 using the "poor man's 
scaling" technique [18], we present exact numerical data of \ f° r specific realizations of the closely 
related Hamiltonian Eq. Q. A key point is that the non-perturbative numerical data enables us to 
quantitatively verify that the perturbative "scaling" approach is sufficient to capture many aspects 
of the mesoscopic fluctuations, as well as to comment on the effect of these fluctuations in the low 
temperature regime where perturbation-theory is not valid. To get a non-perturbative hold on the 
Hamiltonian Eq. ([Q, we use a Quantum Monte-Carlo (QMC) method due to Hirsch and Fye [19]. 
Because only the d-site in the Anderson model is interacting, the introduction of one auxiliary bosonic 
variable renders the entire action quadratic in the fermionic fields. Integrating out the fermions leaves a 
determinant which is a functional of the bosonic variable. Finally, using the determinant as the weight 
to do importance sampling of the bosonic configuration, one can calculate imaginary time correlators 
of the Anderson model. We will study the correlator \ defined in Eq. (|3}. 

Fig-[D a ) shows xCO f° r f° ur different realizations of the conduction electron sea. The first one 
(open symbols) corresponds to the "bulk-like" case, by which we mean a discrete spectrum but without 
mesoscopic fluctuations (constant spacing of the energy levels, no fluctuations of the wave functions). 
The perfect agreement with the universal curve demonstrates that, at the temperature we consider, the 
discreetness of the spectrum per se plays no role^ 1 ) For the three remaining curves (filled symbols), 



( 1 )This implies the absence of any even-odd effects such as the ones observed in [9], as well as any difference between the 
grand canonical and canonical ensembles. At lower T such effects would appear and could be captured by QMC techniques. 
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the one-body matrix describing the conduction electron Hamiltonian is drawn from the Gaussian or- 
thogonal ensemble (GOE), which in particular implies independent Porter- Thomas distributions for 
the intensities |0 Q (O)| 2 . Implicit here is the assumption that — > oo, which, as discussed below, is 
a well-defined limit in the chaotic case. 

An important result arising from these calculations is that not all of the QMC data can be described 
by the universal form, x — fiT/Tx) /T, with a realization dependent Tk- This is particularly evident 
in cases in which \ turns down at low-T since /(T/Tr:) /T is a monotonically decreasing function. A 
more detailed analysis shows that for the parameters used in Fig.^ only about half of the realizations 
yield a temperature dependent susceptibility that can be scaled onto the universal curve with reasonably 
good accuracy. Clearly, mesoscopic fluctuations can cause significant deviations from universality. 

Poor Man 's Scaling. - The basic idea of this technique is that Hamiltonian H with bandwidth 
D cut can be replaced with another Hamiltonian H' with smaller bandwidth while keeping the low 
energy physics unchanged. In the bulk case, the only coupling parameter is Jpo, and its "renormal- 
ization" causes it to diverge as the bandwidth approaches an energy scale that may be identified with 
Tjf . In our case of a discrete spectrum, we take out a single level at a time. Following Anderson, upon 
removing the topmost level cj)p, in order to keep low energy physics unchanged, the amplitude J ai 
between two low energy states changes by <5J Q7 = J a p J^/ep. Interestingly, this equation is solved 
by 5J ai of the form (6 J)<£* (0)^(0) with 

8j=j*\*£mi, ^= j2 ^m. (4) 

ep oD D 

Thus in the "poor man's scaling" approach, the effect of reducing D is captured by letting J flow, 
obtaining at low temperature T the effective coupling constant 

Joff(T) " i jr^ f u, ' (5) 
1 - J J T pi oc (e)de/e 

where J is the "bare" coupling constant defined at the cutoff of the Kondo model Eq. (|2j. In the bulk 
problem it is well known that the "poor man's scaling" result can also be reproduced by two other 
perturbative approaches, namely a one-loop renormalization group analysis, and Abrikosov's method 
for resumming the "parquet" diagrams [20]. This is also true in the presence of fluctuations (see 
e.g. [8] in the perturbative renormalization group context). By working at finite temperature (using 
Matsubara Green functions in conjunction with the other perturbative approaches) |0i oc (e) in Eq. Q is 
replaced by a more natural smoothed LDOS 

/ \ w f°° pioc(e) ,„ 
p sm {uj) = - / de (6) 

n J-oo U + e 

In the following we shall discuss our results in terms of the smoothed LDOS Eq. (|6j. We note that 
one-loop RG analysis has been shown to be a controlled approximation only for a flat (bulk) density 
of states; hence, for the mesoscopic regime, we turn to a comparison with the exact numerics. 

Comparison of Numerics with Poor Man 's Scaling. - At a qualitative level, the relevance of 
p sm (u>) can be assessed by comparing panels (a) and (b) of Fig. ^ where we have plotted the functions 
x(T) and p sm (w) for the same realizations. Note that the cases with large deviations from the universal 
behavior are those for which p sm changes the most from the high to low energy scales. The non- 
monotonic behavior of \ f° r one realization may be understood as being due to a peak in p sm at low 
frequencies. In other cases a dip in p sm causes \ to increase continuously, perhaps until T ~ A. In 
such cases there is no "Fermi-liquid" regime (x ~ constant and T > A) even though the bulk case 
has such a regime. 
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As the Kondo temperature is defined as the scale at which the effective coupling constant J e g 
diverges, it is natural to evaluate the realization dependent fluctuations of Tk around the bulk value 
T K using the perturbative Kondo temperature T£ CIt defined as 

[° Cllt , ^ f^- r s&>> 

/ Psm(w) = / P sm (w) , (7) 

where ~p sm is averaged over realizations and so is by construction essentially constant. The lines in 
Fig-[D a ) correspond to the universal curves x(T) = f(T/T^ crt )/T. Clearly, the high temperature (i.e. 
above T£ crt ) part of the QMC data is well described in this way. To evaluate how quantitatively Eq. 
describes the data, we furthermore consider all our realisations that show approximately "universal" 
behavior, in the sense that the data follow the bulk universal curve obtained from NRG rescaled with 
an appropriate fitted Kondo temperature. This constitutes approximately two-thirds of the realisations 
studied. We have then compared (inset in Fig.0 this TS RG fit with the energy scale 7* crt given by 
Eq. 0. The agreement is quite good, especially considering that Eq. @ is meant to define an energy 
scale rather than an exact number. 

It has to be borne in mind, however, that with the Hirsch and Fye algorithm we are using, the 
ratio U/T that can be practically investigated cannot be made much larger than 500. The parameters 
used in Fig. ^ have been chosen in such a way that temperature significantly smaller than the bulk 
Kondo temperature can be reached, so that the absence of scaling behavior in the low temperature 
regime could be demonstrated. This however means that Z? cu t = 0.18Z7 ~ 16. 6T^- is only a little 
bit more than a decade larger than the bulk Kondo temperature. Here this implies that a good part of 
the fluctuations seen in the high temperature part of Fig.Qcould be explained by fluctuations of the 
smooth density of states p sm at D cut . 

To show that the one-loop renormalization analysis indeed captures the fluctuations associated with 
the energy dependence of the density of states, we have performed another set of QMC simulations, 
with the parameters T = 250A and U = 1000A, so that D cut = 180A and T K ~ 5.1A. The larger 
value of £> cu t already implies that p sm (D cu t) shows significantly less fluctuations. To completely 
exclude the possibility that the observed fluctuation of the susceptibility is an effect of variations 
of the smooth density of states at the scale D cut , we furthermore rescale all the density of states 
[i.e. multiply the (j> a (0) by a constant] so that p sm (-D cu t) = po for all realizations. The resulting 
temperature dependent susceptibility, for four realizations (selected because they display large density 
of states fluctuations), is shown in Fig.|2] The comparison with the universal curves f(T/T^ clt )/T 
shows that in this case, also, the high temperature part of the QMC data is well described in this way. 
It is particularly remarkable, in both Figs. 1 and 2, that the RG result works for temperatures as low as 
Tk or even below. 

We thus reach the central result of this paper: The excellent agreement seen in Fig. between 
theoretical predictions obtained in this way and exact numerical calculations provides a numerical 
demonstration that this approach is a sound starting point for including mesoscopic fluctuations. 

Fluctuations o/Tk- - This agreement justifies using TE 6 as the crossover temperature from 
weak to strong coupling, even for realisations for which the universal scaling is destroyed by the 
mesoscopic fluctuations. Within our RMT modeling, the variance of T^ elt can be easily evaluated 
from Eq. in the limit 8T K = (T^ 01 * — T K ) <C T K as [21] 

(STk 2 ) „ A 

For our model, a ~ (4/ir) In 2; more generally it is a number of order 1. An improved model can 
be developed through semiclassical analysis, in which case (■) may signify an energy average. The 
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Fig. 2 - Triangles: QMC calculation of the local susceptibility for four realizations, with the parameters T = 250 A 
and U = 1000A, so that D cut = 180A ~ 35. 3T£, with T% ~ 5.1A. Lines: universal curves using T£ ert from 
Eq. 0. Inset: Corresponding smoothed densities of states (with matching line styles). 



correlator (p sm (u>)p S m(w')) is then related to the Fourier transform of the classical probability P c \(t) 
that a particle starts at the impurity with energy Ep and returns there after time t [22]. Eq. (|8} can 
therefore be generalized to 



(9) 



- K 

(V is the volume of the dot). For chaotic systems, 



P c i0) = for t < h/E T h (10) 

P c i(i) = l/V for t > h/E Th 

so that Eq. l|8} is then trivially recovered in the limit y = E^/T^ — > oo. The fluctuations of Tk 
vanish, of course, in the opposite limit y — > 0. A remarkable implication of Eqs. d9l>- (llOI > is that the 
value of ETh does not affect the mesoscopic fluctuations for a chaotic system with i?Th ^> Tk. The 
fluctuations are in this sense "universal", and are described by RMT for any chaotic system. 

Conclusions. - We have shown that the mesoscopic fluctuations typical in nanostructures signif- 
icantly affect the physics of a Kondo impurity. Our main result is a demonstration that the fluctuations 
of the high temperature properties can be understood by following a simple poor man's scaling argu- 
ment (or, equivalently, one-loop RG or resummation of parquet diagrams) and that this works down to 
T w Tk- In addition, we show (1) deviations from universal behavior at low temperature, (2) an ex- 
plicit formula for the variance of the Kondo temperature, and (3) a semiclassical connection between 
the fluctuations of Tk and the classical probability of return. 

Kondo physics has recently received renewed attention because of its relevance to the low tem- 
perature physics of quantum dots [14, 15]. These "artificial magnetic impurities" provide tunability 
of the microscopic Hamiltonian, and can be naturally interfaced with larger quantum dots displaying 
the kind of mesoscopic fluctuations that we address here. From an experimental point of view, the 
physics we are discussing should be easily observable in such systems. In particular, we stress that 
though we have here specifically considered the local susceptibility, the qualitative features we have 
discussed - the mere existence of the mesoscopic fluctuations, the absence of universality of the low 
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temperature regime, and the universal character, but with a realisation dependent Kondo temperature, 
of the T > Tk regime - should be observed for other physical quantities. In the same way, the fact 
that the the realisation dependent Kondo temperature can be extracted from a fit to the universal form 
in the high temperature regime implies that our quantitative predictions Eqs. (|8j-(|9j can also be tested 
on physical quantities other than the local susceptibility. 

Finally, we mention that experimental set-ups where the correlator Eq. (0 can be directly measured 
are also possible. In particular, the problem of charge fluctuations in a gated, Coulomb blockaded 
metal grain [23] has been shown to be equivalent to a Kondo problem [24]. In this mapping the charge 
and gate voltage of the CB problem map to the spin of the Kondo impurity and a local magnetic field 
applied to it. The differential capacitance at a degeneracy point in the CB problem is thus equivalent 
to the local susceptibility Eq. of a Kondo problem. In the presence of a strong polarizing magnetic 
field, the equivalent Kondo Problem is described by exactly our Hamiltonian Eq. (0. Thus, a gated dot 
connected to another dot provides a realization of a Kondo impurity coupled to a mesoscopic electron 
gas where it is possible to measure the correlator Eq. directly. 

* * * 
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